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The event-by-event fluctuations in heavy ion collisions carry information about the 
thermodynamic properties of the hadronic system at the time of freeze-out. By 
studying these fluctuations as a function of varying control parameters, such as the 
collision energy, it is possible to learn much about the phase diagram of QCD. As a 
timely example, we stress the methods by which present experiments at the CERN 
SPS can locate the second order critical point at which a line of first order phase 
transitions ends. Those event-by-event signatures which are characteristic of freeze- 
out in the vicinity of the critical point will exhibit nonmonotonic dependence on 
control parameters. We focus on observables constructed from the multiplicity and 
transverse momenta of charged pions. We find good agreement between NA49 data 
and thermodynamic predictions for the noncritical fluctuations of such observables. 
We then analyze the effects due to the critical fluctuations of the sigma field. We 
estimate the size of these nonmonotonic effects which appear near the critical point, 
including restrictions imposed by finite size and finite time. 



1 Introduction 

In my talk at Strong and Electroweak Matter '98, I presented recent work 
on the physics which arises in two different areas of the QCD phase diagram. 
Cold dense quark matter forms a color superconductor, and I compared the 
superconducting phase expected in QCD with two massless quarks with that 
expected for three massless quarks in which chiral symmetry is broken by 
color-flavor locking. Alford, Berges and I have recently completed an analysis 
of the phase diagram of zero temperature QCD as a function of density and 
strange quark mass, in Ref. □. This brings the ideas I presented in Copenhagen 
together into a consistent picture, and I refer you to that paper for an up-to- 
date treatment of the subject, and references to the literature. 

The other half of my talk in Copenhagen was a sketch of methods by which 
present heavy ion experiments can find the critical point in the QCD phase 
diagram at nonzero temperature T and baryon chemical potential /i. Like the 
end point in the electroweak phase diagram, discussed by others at this meet- 
ing, this critical point is a second order transition in the Ising universality class 
which occurs at the end of a line of first order phase transitions. Stephanov, 

"Work done in collaboration with Misha Stephanov and Edward ShuryakE] 
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Shuryak and I have recently completed a detailed analysis of the signatures 
of the physics characteristic of the vicinity of this pointl^ begun in Ref. u. I 
described this work in a preliminary fashion in Copenhagen; in these proceed- 
ings, I summarize the results and implications, of Ref. p. Those interested in 
the derivation of these results should see Ref. □. 

Large acceptance detectors, such as NA49 and WA98 at CERN, have made 
it possible to measure important average quantities in single heavy ion collision 
events. For example, instead of analyzing the distribution of charged particle 
transverse momenta obtained by averaging over particles from many events, we 
can now study the event-by-event variation of the mean transverse momentum 
of the charged pions in a single event, Although much of this data still 

has preliminary status, with more statistics and more detailed analysis yet to 
come, some general features have already been demonstrated. In particular, 
the event-by-event distributions of these observables are as perfect Gaussians 
as the data statistics allow, and the fluctuations — the width of the Gaussians 
— are smallQ 

This is very different from what one observes in pp collisions, in which 
fluctuations are large. These large non-Gaussian fluctuations clearly reflect 
non-trivial quantum fluctuations, all the way from the nucleon wave function 
to that of the secondary_hadrons, and are not yet sufficiently well understood. 
As discussed in Refs. Era, thermal equilibration in AA collisions drives the 
variance of the event-by-event fluctuations down, close to the value determined 
by the variance of the inclusive one-particle distribution divided by the square 
root of the multiplicity. 

Can we learn something from the magnitude of these small fluctuations 
and their dependence on the parameters of the collision? What do the widths 
of the Gaussians tell us about the thermodynamics of QCD? Some of these 
questions have been addressed in Refs. 019 where it was pointed out that, for 
example, temperature fluctuations are related to heat capacity via 

((AT) 2 ) _ 1 



T 2 C V (T)> 



(1) 



and so can tell us about thermodynamic properties of the matter at freeze-out. 
Furthermore, Mrowczyhski has discussed the study of the compressibility of 
hadronic. matter at freeze=put via the event-bj^event fluctuations of the particle 
numberta and GazdzickiL3 and Mrowczyhski liil have considered event-by-event 
fluctuations of the kaon to pion ratio as measured by NA49Q. In pp physics 
one can hope to extract quantum mechanical information about the initial 

b We denote the mean transverse momentum of all the pions in a single event by px rather 
than (pt) because we choose to reserve (. . .} for averaging over an ensemble of events. 
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state from event-by-event fluctuations of the final state; in heavy ion collisions 
equilibration renders this an impossible goal. In AA collisons, then, the new 
goal is to use the much smaller, Gaussian event-by-event fluctuations of the 
final state to learn about thermodynamic properties at freeze-out. 

It is worth noting that once a large acceptance detector has presented 
convincing evidence that the event-by-event distribution of, for example, pt 
is Gaussian, then the measurement of the width of such a distribution can 
be accomplished by "event-by-event" measurements in which only two pions 
per escent are observed. This has recently been emphasized by Bialas and 
KochtL3 Of course, this approach measures the width of the event-by-event 
distribution whether or not it is Gaussian; it is only the results of a large 
acceptance experiment like NA49 which motivate a thermodynamic analysis 
of the event-by-event fluctuations. 

Stephanov, Shuryak and I focus on observables constructed from the multi- 
plicity and the momenta of the charged particles in the final state, as measured 
by NA49. We leave the extension of the methods of this paper to the study of 
thermodynamic implications of the NA49 Gaussian distribution of event-by- 
event K/tt ratios □ and of the WA98 Gaussian distribution of event-by-event 
7T° /n^ ratiosE^I for future work. 

One of the lessons of our paper is that it is difficult to apply thermodynamic 
relations like |l]) directly. To see a sign of this, note that the event-by-event 
fluctuations of the energy E of a part of a finite system in thermal equilibrium 
are given by ((A_E) 2 ) = T 2 Cy(T). For a system in equilibrium, the mean 
values of T and E are directly related by an equation of state E(T); their 
fluctuations, however, have quite different behavior as a function of Cy, and 
therefore behave differently when Cy diverges at a critical point. The fluctua- 
tions of "mechanical" observables increase at the critical point. Because T(E) 
is singular at the critical point, the fluctuations of T decrease there. It is a fact 
that what we measure are the mechanical observables, and since we in general 
only know T{E) for simple systems we call thermometers, we cannot apply (0) 
to the complicated system of interest. It is not in fact necessary to translate 
the observed "mechanical" variable (the mean transverse momentum pt for 
example) into a temperature in order to detect the critical point. It is easier 
to look directly at the fluctuations of observable quantities. We demonstrate 
that the fluctuations of px grow at the critical point. 

Although our methods are general, we focus in Ref. on how to use them 
to find and study the critical end-point E on the phase diagram of QCD in 
the T/z plane. The possible existence of such a point, as an endpoint of the 
first order transition separating quark-gluon plasma from hadron matter, and 
its universal critical properties have been pointed out recently in Refs. E-aO. 
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The point E can be thought of as a descendant of a tricritical point in the 
phase diagram for 2-flavor QCD with massless quarks. In a previous, letter, 
we have laid out the basic ideas for observing the critical endpoint u. The 
signatures proposed in Ref. are based on the fact that such a point is a 
genuine thermodynamic singularity at which susceptibilities diverge and the 
order parameter fluctuates on long wavelengths. The resulting signatures all 
share one common property: they are nonmonotonic as a function of an ex- 
perimentally varied parameter such as the collision energy, centrality, rapidity 
or ion size. Once experimentalists vary a control parameter which causes the 
freeze-out point in the (T, /u,) plane to move toward, through, and then past the 
vicinity of the endpoint E, they should see all the signatures we describe first 
strengthen, reach a maximum, and then decrease, as a nonmonotonic function 
of the control parameter. It is important to have a control parameter whose 
variation changes the fi at which the system crosses the transition region and 
freezes out. The collision energy is an obvious choice, since it is known experi- 
mentally that varying the collision energy has a large effect on /i at freeze-out. 
Other possibilities should also be explored.^ 

We assume throughout that freeze-out occurs from an equilibrated hadronic 
system. If freeze-out occurs "to the left" (lower /j,; higher collision energy) of 
the critical end point E, it occurs after the matter has traversed the crossover 
region in the phase diagram. If it occurs "to the right" of E, it occurs after the 
matter has traversed the first order phase transition. This is the situation in 
which our assumption of freeze-out from an equilibrated system is most open to 
question. First, one may imagine hadronization directly from the mixed phase, 
without time for the hadrons to rescatter. Hadronic elastic scattering cross- 
sections are large enough that this is unlikely. Second, one may worry that the 
matter is inhomogeneous after the first order transition, and has not had time 
to re-equilibrate. Fortunately, our assumption is testable. If the matter were 
inhomogeneous atfreeze-out, one can expect non-Gaussian fluctuations in var- 
ious observablesliS which would be seen in the same experiments that seek the 
signatures we describe. We focus on the Gaussian thermal fluctuations of an 
equilibrated system, and study the nonmonotonic changes in these fluctuations 
associated with moving the freeze-out point toward and then past the critical 
point, for example from left to right as the collision energy is reduced. 

Ref. is devoted to a detailed analysis of the physics behind event-by-event 
fluctuations in relativistic heavy ion collisions and the resulting effects unique 
to the vicinity of the critical point in the phase diagram of QCD. Most of 

c If the system crosses the transition region near E, but only freezes out at a much lower 
temperature, the event-by-event fluctuations will not reflect the thermodynamics near E. In 
this case, one can push freeze-out to earlier times and thus closer to E by using smaller ionsu 
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our analysis is applied to the fluctuations of the observables characterizing the 
multiplicity and momenta of the charged pions in the final state of a heavy ion 
collision. There are several reasons why the pion observables are most sensitive 
to the critical fluctuations. First, the pions are the most numerous hadrons 
produced and observed in relativistic heavy ion collisions. A second, very 
important reason, is that pions couple strongly to the fluctuations of the sigma 
field (the magnitude of the chiral condensate) which is the order parameter of 
the phase transition. Indeed, the pions are the quantized oscillations of the 
phase of the chiral condensate and so it is not surprising that at the critical end 
point, where the magnitude of the condensate is fluctuating wildly, signatures 
are imprinted on the pions. 

2 Noncritical Thermal Fluctuations in Heavy Ion Collisions 

Before we discuss the effects of the critical fluctuations, we must analyze the 
thermal fluctuations which are present if freeze-out does not occur in the vicin- 
ity of the critical point. In this section, but not throughout this paper, we 
assume that the system freezes out far from the critical point in the phase 
diagram, and can be approximated as an ideal resonance gas when it freezes 
out. We compare some of our results to preliminary data from the NA49 ex- 
periment on PbPb collisions at 160 AGeV, and find broad agreement. The 
results obtained seem to support the hypotheses that most of the fluctuation 
observed in the data is indeed thermodynamic in origin, and that this system 
is not freezing out near the critical point. 

As a first test of our resonance gas model, we analyze the fluctuations in 
an ideal Bose gas of pions, and then add as many of the effects which this 
simple treatment neglects except that we assume that no effects due to critical 
fluctuations are significant. We model the matter in a relativistic heavy ion 
collision at freeze-out as a resonance gas in thermal equilibrium, and begin by 
calculating the variance of the event-by-event fluctuations of total multiplicity 
N. The fluctuations in N are not affected by the boost which the pion momenta 
receive from the collective flow, but they are contaminated experimentally 
by fluctuations in the impact parameter. This experimental contamination 
can be reduced by making a tight enough centrality cut using a zero degree 
calorimeter. 

We find {(AN) 2 )/(N) w 1.5, which we compare with NA49 results from 
central Pb-Pb collisions at 160 AGeV. It is clear that with no cut on centrality, 
one would see a very wide non-Gaussian distribution of multiplicity determined 
by the geometric probability of different impact parameters b. Gaussian ther- 
modynamic fluctuations can only be seen if a tight enough cut in centrality 
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is applied. The event-by-event iV-distribution found by NA49 when they use 
only the 5% most central of all events, with centrality measured using a zero 
degree calorimeter, is Gaussian to within about 5%. This cut corresponds to 
keeping collisions with impact parameters b < 3.5 fmQ The non-Gaussianity 
could be further reduced by tightening the centrality cut further. From the 
data, we have ((AN) 2 )/(N) = 2.008 ± 0.009, which suggests that about 75% 
of the observed fluctuation is thermodynamic in origin. The contamination 
introduced into the data by fluctuations in centrality could be reduced by an- 
alyzing data samples with more or less restrictive cuts but the same (N), and 
extrapolating to a limit in which the cut is extremely restrictive. This could be 
done using cuts centered at any centrality. Our resonance gas model predicts 
that as the centrality cut is tightened, the ratio v 2 hc (N)/(N) should decrease 
toward a limit near 1.5. 

Although further work is certainly required, it is already apparent that the 
bulk of the multiplicity fluctuations observed in the data are thermodynamic 
in origin, impact parameter. Note that our prediction is strongly dependent 
on the presence of the resonances; had we not included them, our prediction 
would have been significantly lower, farther below the data. Because the multi- 
plicity fluctuations are sensitive to impact parameter fluctuations, it may prove 
difficult to explain their magnitude with greater precision even in future. How- 
ever, the fact that they are largely thermodynamic in origin suggests that the 
effects present near the critical point, which we describe below, could result in 
a significant nonmonotonic enhancement of the multiplicity fluctuations. This 
would be of interest whether or not the noncritical fluctuations on top of which 
the nonmonotonic variation occurs are understood with precision. 

We then turn to a calculation of the variance of the event-by-event fluctu- 
ations of the mean transverse momentum, px- We first calculate the width of 
the inclusive pT-distribution, Vinc(pT). In the absence of any correlations, the 
event-by-event fluctuations of the mean transverse momentum of the charged 
pions in an event, Vebcipr) = ((Apr) 2 ) 1 / 2 would be given by Vi nc (pT)/ {N) 1 / 2 , 
and this turns out to be a very good approximation in the present data as 
we discuss below. We calculate numerically the contribution to t>i nc (PT) from 
"direct pions", already present at freeze-out, and from the pions generated 
later by resonance decay. We have simulated a gas of pions, nucleons and res- 
onances in thermal equilibrium at freeze-out, including the 7r, K, r/, p, u>, rj', 
N, A, A, £ and S, and then simulated the subsequent decay of the resonances. 
That is, we have generated an ensemble of pions in three steps: (i) Thermal 
ratios of hadron multiplicities .were calculated assuming equilibrium ratios at 
chemical freeze-out. FollowingO, the values T c h = 170 MeV and /ibaryon = 200 
MeV were used, (ii) Then, a program generates hadrons with multiplicities 
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determined at chemical freezeout, but with thermal momenta as appropriate 
at the thermal freeze-out temperature, which we take to be Tf = 120 MeV, 
with [in = 60 MeV. The last step (iii) is to decay all the resonances. From 
the resulting ensemble of pions (the sum of the direct pions and those from 
the resonances) we obtain Vi nc (pT) / (pt) = 0.66 The resonances turn out to be 
less important here than in the calculation of the multiplicity fluctuations, in 
that the resonance gas prediction for v- lnc {pT) / '(pt) is almost indistinguishable 
from that of an ideal Bose gas of pions. 

To this point, we have calculated the fluctuations in px as if the matter 
in a heavy ion collision were at rest at freeze-out. This is not the case: by 
that stage the hadronic matter is undergoing a collective hydrodynamic ex- 
pansion in the transverse direction, and this must be taken into account in 
order to compare our results with the data. A very important point here is 
that the fluctuations in pion multiplicity are not affected by flow, and our 
prediction for them is therefore unmodified. However the event-by-event fluc- 
tuations of mean px are certainly affected by flow. The fluctuations we have 
calculated pertain to the rest frame of the matter at freeze-out, and we must 
now boost them. A detailed account of the resulting effects would require 
a complicated analysis. We use the simple approximation Ii3 that the effects 
of flow on the pion momenta can be treated as a Doppler blue shift of the 
spectrum: ti{pt) — * ti{pt\/^ — + ft)- This blue shift increases (pt), and 

increases v- mc {pT), but leaves the ratio Vi nc (pT) / (pt) (and therefore the ratio 
Vebeipr) / (pt)) unaffected. However, event-by-event fluctuations in the flow 
velocity p must .still be taken into account. This issue was discussed qualita- 
tively already inn, where it was argued that this effect must be relatively weak. 
In Ref. Id we provide the first rough estimate of its magnitude. We estimate 
that fluctuations in the flow velocity increase Vi nc (j>T) / ' {pt) from 0.66 to 0.67. 
The largest uncertainty in our estimate for Uin C (fr) / ' (jpt) is not due to the 
fluctuations in the flow velocity, which can clearly be neglected, but is due to 
the velocity itself. The blue shift approximation which we have used_applies 
quantitatively only to pions with momenta greater than their mass El Be- 
cause of the nonzero pion mass, boosting the pions does not actually scale the 
momentum spectrum by a momentum independent factor. Furthermore, in a 
real heavy ion collision there will be a position dependent profile of velocities, 
rather than a single velocity f3. A more complete calculation of Uin C (pr) / '(pr) 
would require a better treatment of these effects in a hydrodynamic model; we 
leave this for the future. 

We compare our results to the NA49 data, in which v- mc (pT) / '(pt) = 
0.749 ± 0.001. We see that the major part of the observed fluctuation in pt 
is accounted for by the thermodynamic fluctuations we have considered. Our 
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prediction is about 10% lower than that in the data. First, this suggests that 
there may be a small nonthermodynamic contribution to the pT-fluctuations, 
for example from fluctuations in the impact parameter. (However, we expect 
that the fluctuations of an intensive quantity like pr are less sensitive to im- 
pact parameter fluctuations than are those of the multiplicity, and this seems 
to be borne out by the data.) The other source of the discrepancy is the blue 
shift approximation. We leave a more sophisticated treatment of the effects of 
flow on the spectrum to future work. Such a treatment is necessary before we 
can estimate how much of the 10% discrepancy is introduced by the blue shift 
approximation. Future work on the experimental side (varying the centrality 
cut) could lead to an estimate of how much of the discrepancy is due to impact 
parameter fluctuations. 

We have gone as far as we will go in this paper in our quest to understand 
the thermodynamic origins of the width of the inclusive single particle distri- 
bution. We now turn to the ratio of the scaled event-by-event variation to the 
variance of the inclusive distribution: 

VF=— , ' = 1.002 ±0.002. 2 

finc(PT) 

The difference between the scaled event- by- event variance and the variance of 
the inclusive distribution is less than a percent in the NA49 data.^ 

We analyze a number of noncritical contributions to the ratio V~F, which 
we write 

Vf = Vf b f ics f ec . (3) 

Fb is the contribution of the Bose enhancement of the fluctuations of identical 
pions. We calculate this effect and find y/Fs ~ 1-02. F ICS describes the effect 
of the correlations induced by the fact that pions produced by the decay of 
a resonance after freeze-out do not have a chance to rescatter. We estimate 
it by dividing the pions from our resonance gas simulation into "events" of 
varying sizes, and evaluating F. Since Bose enhancement is not included in the 

d We explain in an Appendix in Ref. Q that in order to be sure that F = 1 when there 
are no correlations between pions, care must be taken in constructing an estimator for 
v„\ h , (pt) using a finite sample of events, each of which has finite multiplicity. The appropriate 
prescription is to weight events in the event-by-event average by thejj multiplicity, and we 
have made the appropriate correction in writing (H). Other authors have introduced the 
correlation measure <t> PT = (N) 1 / 2 v e ^ te (pj') — t>i n c(pT)- Because Ui n c(PT) i s scaled by the 
blue shift introduced by the expansion velocity, so is <3> PT . This makes $ PT harder to 
predict than F. However, for convenience, we note that if one uses the experimental value 
of «i nc (pT)i a value y/F = 1.01 corresponds to <E> PT = 2.82 MeV, and the \fF in the data 
(|) corresponds to <S> PT = 0.6 ± 0.6 MeV. 
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simulation, the F so obtained is just F ros . We find no statistically significant 
contribution, and conclude that |F rcs — 1| < 0.01. 

The third contribution, -FeCj is due to energy conservation in a finite sys- 
tem. This is most easily described by considering the event-by-event fluctu- 
ations An p in the number of pions in a bin in momentum space centered at 
momentum p. Consider the correlator (An p Arik) ■ When one n p fluctuates 
up, others must fluctuate down, and it is therefore more likely that fluctu- 
ates downward. Energy conservation in a finite system therefore leads to an 
anti-correlation which is off-diagonal in pk space. u e b e (f>r) is determined by 
(An p Ank), and the result of this anti-correlation is a reduction: 

V^-0.99. (4) 

If the observed charged pions are in thermal contact with an unobserved heat 
bath, the anti-correlation introduced by energy conservation decreases as the 
heat capacity of the heat bath increases. The estimate (||) assumes that the 
heat capacity of the direct charged pions is about 1/4 of the total heat ca- 
pacity of the hadronic system at freezeout. In addition to the contributions 
we calculate, is affected by the finite two-track resolution in the detector, 
and by final state Coulomb interactions between charged pions. NA49 esti- 
mates that these contributions reduce F by about the same amount that Bose 
enhancement increases it. 

We conclude that the ratio v~F measured by NA49 is broadly consistent 
with thermodynamic expectations. It receives a positive contribution from 
Bose enhancement, negative contributions from energy conservation and two- 
track resolution, and a positive contribution from the effect of resonance de- 
cays. These contributions to VF are all roughly at the 1% level (or smaller 
in the case of that from resonance decays) and it seems that they cancel in 
the data (0). Our results support the general idea that the small fluctua- 
tions observed in ^4^4 collisions, relative to those in pp, are consistent with 
the hypothesis that the matter in the AA collisions achieves approximate local 
thermal equilibrium in the form of a resonance gas. 

With more detailed experimental study, either now at the SPS, or soon at 
RHIC (STAR will study event-by-event fluctuations in px, N, particle ratios, 
etc; PHENIX and PHOBOS in N only) it should be possible to disentangle the 
different effects we describe. Making a cut to look at only low px pions should 
increase the effects of Bose enhancement. The anti-correlation introduced by 
energy conservation is due to terms in (An p Ank) which are off-diagonal in 
pk. Thus, a direct measurement of (An p Ank) would make it easy to separate 
this anti-correlation from other effects. The cross correlation (AN Apr) is 
also a very interesting observable to study. It vanishes for a classical ideal 
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gas. This means that whereas v c ^ c (pt) receives a dominant contribution from 
the width of the inclusive single particle distribution, this effect cancels in 
(ANApt) and the remaining effects due to Bose enhancement and energy 
conservation dominate. Although this cross-correlation is small, it is worth 
measuring because it only receives contributions from interesting effects. 

We hope that the combination of the theoretical tools we have provided 
and the present NA49 data provide a solid foundation for the future study of 
the thermodynamics of the hadronic matter present at freeze-out in heavy ion 
collisions. Once data is available for other collision energies, centralities or ion 
sizes, the present NA49 data and the calculations of this section will provide an 
experimental and a theoretical baseline for the study of variation as a function 
of control parameters. 

Our analysis demonstrates that the observed fluctuations are broadly con- 
sistent with thermodynamic expectations, and therefore raises the possibility 
of large effects when control parameters are changed in such a way that ther- 
modynamic properties are changed significantly, as at a critical point. The 
smallness of the statistical errors in the data also highlights the possibility 
that many of the interesting systematic effects we analyze in this paper will be 
accessible to detailed study as control parameters are varied. 

3 Pions Near the Critical Point: Interaction with the Sigma Field 

With the foundations established, we now describe how the fluctuations we 
analyze will change if control parameters are varied in such a way that the 
baryon chemical potential at freeze-out, /if, moves toward and then past the 
critical point in the QCD phase diagram at which a line of first order transitions 
ends at a second order endpoint. The good agreement between the noncritical 
thermodynamic fluctuations we analyze in Section 2 and NA49 data make it 
unlikely that central PbPb collisions at 160 AGeV freeze out near the critical 
point. Estimates we have made in Ref. a suggest that the critical point is 
located at a baryon chemical potential /i such that it will be found at an 
energy between 160 AGeV and AGS energies. This makes it a prime target for 
detailed study at the CERN SPS by comparing data taken at 40 AGeV, 160 
AGeV, and in between. If the critical point is located at such a low /j, that the 
maximum SPS energy is insufficient to reach it, it would then be in a regime 
accessible to study by the RHIC experiments. We want to stress that we are 
more confident in our ability to describe the properties of the critical point 
and thus to predict how to find it than we are in our ability to predict where 
it is. 

We now describe how the fluctuations of the pions will be affected if the 
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system freezes out near the critical endpoint. First, because the pions at freeze- 
out are now in contact with a heat bath whose heat capacity diverges at the 
critical point, the effects of energy conservation parametrized by Fec — 1 arc 
greatly reduced. However, since Fec is close to one even away from the critical 
point, this is a small effect. 

The dominant effects of the critical fluctuations on the pions are the direct 
effects occuring via the ennr coupling. In the previous section, we made the 
assumption that the "direct pions" at freeze-out could be described as an ideal 
Bose gas. We do not expect this to be a good approximation if the freeze-out 
point is near the critical point. The sigma field is the order parameter for the 
transition and near the critical point it therefore develops large critical long 
wavelength fluctuations. These fluctuations are responsible for singularities in 
thermodynamic quantities. We find that because of the Gamr coupling, the 
fluctuations of both the multiplicity and the mean transverse momentum of 
the charged pions do in fact diverge at the critical point. 

We then estimate the size of the effects in a heavy ion collision. This 
requires first estimating the strength of the coupling constant G, and then 
taking into account the finite size of the system and the finite time during which 
the long wavelength fluctuations can develop. We find a large increase in the 
fluctuations of both the multiplicity and the mean transverse momentum of the 
pions. This increase would be divergent in the infinite volume limit precisely 
at the critical point. We apply finite size and finite time scaling to estimate 
how close the system created in a heavy ion collision can come to the critical 
singularity, and consequently how large an effect can be seen in the event-by- 
event fluctuations of the pions. We conclude that the nonmonotonic changes 
in the variance of the event-by-event fluctuation of the pion multiplicity and 
momenta which are induced by the universal physics characterizing the critical 
point can easily be between one and two orders of magnitude greater than the 
statistical errors in the present data. 

The value of the coupling G in vacuum can be estimated either from the 
relationship between the sigma and pion masses and f v or from the width of the 
sigma. Both yield an estimate G ~ 1900 MeV, where we have used nig- = 600 
MeV. The width of the sigma is so large that this "particle" is only seen as 
a broad bump in the s-wave ir — tt scattering cross-section. The vacuum ennr 
coupling must be at least as large as G <~ 1900 MeV, since the sigma would 
otherwise be too narrow. 

The vacuum value of G would not change much if one were to take the 
chiral limit m — ► 0. The situation is different at the critical point. Taking 
the quark mass to zero while following the critical endpoint leads one to the 
tricritical point P in the phase diagram for QCD with two massless quarks. At 
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this point, G vanishes as we discuss-below. This suggests that at E, the coupling 
G is less than in vacuum. In Ref. u, we use what we know about physics near 
the tricritical point P to make an estimate of how much the coupling G is 
reduced at the critical endpoint E (with the quark mass m having its physical 
value), relative to the vacuum value G ~ 1900 MeV estimated above. ln 

We begin by recalling some known results. (For details, see Refs. tSllafi.) 
In QCD with two massless quarks, a spontaneously broken chiral symmetry 
is restored at finite temperature. This transition is likely second order and 
belongs in the universality class of 0(4) magnets in three dimensions. At zero 
T, various models suggest that the chiral symmetry restoration transition at 
finite /i is first order. Assuming that this is the case, one can easily argue that 
there must be a tricritical point P in the Tfi phase diagram, where the transi- 
tion changes from first order (at higher than P) to second order (at lnwpr ^), 
and such a tricritical point has been found in a variety of modelsEJ'Ejnj The 
nature of this point can be understood by considering the Landau-Ginzburg 
effective potential for (j> a , order parameter of chiral symmetry breaking: 

fK0a) = -j$a$a + -^{4>a4>a) 2 + g {4>a4>a ) 3 ■ (5) 

The coefficients a, b and c > are functions of [i and T. The second order phase 
transition line described by a = at b > becomes first order when b changes 
sign, and the tricritical point P is therefore the point at which a — h = 0. The 
critical properties of this point can be inferred from universality Oc3, and the 
exponents are as in the mean field theory (||). We will use this below. Most 
important in the present context is the fact that because (4>) = at P, there 
is no CT7T7T coupling, and G = there. 

In real QCD with nonzero quark masses, the second order phase transition 
becomes a smooth crossover and the tricritical point P becomes E, the second 
order critical endpoint of a first order phase transition line. Whereas at P there 
are four massless scalar fields undergoing critical long wavelength fluctuations, 
the <7 is the only field which becomes masslesiL.at the point E, and the point 
E is therefore in the Ising universality classlij'E£l. The pions remain massive at 
E because of the explicit chiral symmetry breaking introduced by the quark 
mass m. Thus, when we discuss physics near E as a function of /i and T, 
but at fixed m, we will use universal scaling relations with exponents from the 
three dimensional Ising model. Our present purpose, however, is to imagine 
varying m while changing T and fi in such a way as to stay at the critical point 
E, and ask how large G (and m w ) become once m is increased from zero (the 
tricritical point P at which G — = 0) to its physical value. For this task, we 
use exponents describing universal physics near P. Applying tricritical scaling 
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relations all the way up to a quark mass which is large enough that m v is not 
small compared to T c may introduce some uncertainty into our estimate. 

We first determine the trajectory of the critical line of Ising critical points 
E as a function of quark mass mj^| and then find that G ~ m 3 / 5 along this line, 
where m is the light quark mass. Thus the coupling G is suppressed compared 
to its "natural" vacuum value G vac by a factor of order (raj Aqcd) 3 ^ 5 • Taking 
Aqcd ~ 200 MeV, m ~ 10 MeV we obtain our estimate 

G E ~ -Jp ~ 300 MeV . (6) 

The main source of uncertainty in this estimate is our inability to compute the 
various nonuniversal masses which enter the estimate as prefactors in front of 
the m dependence which we have followed. In other words, we do not know 
the correct value to use for Aqcd in the suppression factor which we write as 
(WA QCD ) 3/5 . 

The final ingredient we need is an estimate of the correlation length £ of 
the sigma field, which is infinite at the critical point. In practice, there are 
important restrictions on how large £ can become. Two particle interferomctry 
lia suggests that the size of regions over which freeze-out is homogeneous is 
roughly 12 fm in both the longitudinal and transverse directions. This means 
that the finite size of the system limits £ to be less than about this value. The 
finite time restriction is stricter, but harder to estimates Although the size 
of the system allows the correlation length to grow to 12 fm, there may not 
be enough time for such long correlations to grow. We use £ max ~ 6 fm as a 
rough estimate of the largest correlation length possible if control parameters 
are chosen in such a way that the system freezes out close to the critical point. 

We now return to our discussion of the effects of the long wavelength 
sigma fluctuations on the fluctuations of the pions. We use mean field theory 
throughout Ref. p. The fluctuations of the sigma field around the minimum 
of {1(a) are not small; however, this does not make much difference to the 
quantities of interest, all of which diverge like m^ 2 ~ £ 2 at the critical point. 
The divergence is that of the sigma field susceptibility, and for the 3d-Ising 
universality class we know the corresponding exponent to be j/u = 2 — rj 
which is rs 2 to within a few percent because rj is small. We can therefore 
safely use mean-field mean field results with their m~ 2 divergence, and will 
take Too- ~ l/£max ~ 1/(6 fm) in our estimates. 

We now have all the ingredients in place to present our estimate of the size 
of the effect of the critical fluctuations of the sigma field on the fluctuations 

e See Ref. for a derivation of the analogous line of Ising points emerging from the tricritical 
point in the QCD phase diagram at zero /iasa function of m and the strange quarij mass 
m 3 . This tricritical point can be related to the one we are discussing by varying m s B 
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of the direct pions, via the coupling G. We express the size of the effect of 
interest by rewriting the ratio vF of (Q) and (||) as 

VF = ^F B F les F EC F a 

and presenting F a . We find: 

t-, 1 , r, or / ^freeze — out \ / ^freeze— out \ f —\ 

F - = 1 + °- 35 U00MeV ) l^fin-J ' (?) 

where we have taken T — 120 MeV and \i v = 60 MeV. F a will be reduced 
by about a factor of two, because not all of the pions which are observed are 
direct. The coupling G transmits the effects of the critical a fluctuations to 
the pions at freezeout, not to the (heavier) resonances. The size of the effect 
depends quadratically on the coupling G. We argued above that G is reduced 
to Ge ~ 300 MeV at the critical point. However, freeze-out may occur away 
from the critical point, in which case G would be larger, although still much 
smaller than its vacuum value. The size of the effect also depends quadratically 
on the sigma correlation length at freeze-out, and we have seen that there are 
many caveats in an estimate like ^freeze-out ~ £max ~ 6 fm. 

We have studied two different effects of the critical fluctuations on y/F. 
First, .Fec - * 1) leading to about a 1% increase in y/F. The direct effect of 
the critical fluctuations is a much larger increase in y/~F by a factor of y/F^. 
We have displayed the various uncertainties in the factors contributing to our 
estimate (0) so that when an experimental detection of an increase and then 
subsequent decrease in y/~F occurs, as control parameters are varied and the 
critical point is approached and then passed, we will be able to use the mea- 
sured magnitude of this nonmonotonic effect to constrain these uncertainties. 
It should already be clear that an effect as large as 10% in \[F\, is easily pos- 
sible; this would be 50 times larger than the statistical error in the present 
data. 

We now give a brief account of the effect of critical fluctuations on ((AiV) 2 ) 
and (ANApx)- The contribution of the direct pions to ((AN) 2 ) can easily 
double, but the multiplicity fluctuations are dominated by the pions from res- 
onance decay, so we estimate that the critical multiplicity fluctuations lead to 
about a 10-20% increase in ((AN) 2 ). (This neglects the pions from sigma de- 
cay. See below.) The cross-correlation (AN Apr) only receives contributions 
from nontrivial effects, and we find that near the critical point, the contribu- 
tion from the interaction with the sigma field is dominant. We estimate that 
(for Gfreozo-out ~ 300 MeV and ^freeze— out ~ 6 fm) the cross-correlation will be 
a factor of 10-15 times larger than in the absence of critical fluctuations! The 
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lesson is clear: although this correlation is small, it may increase in magnitude 
by a very large factor near the critical point. 

The effects of the critical fluctuations can be detected in a number of ways. 
First, one can find a nonmonotonic increase in the suitably normalized 
increase in the variance of event-by-event fluctuations of the mean transverse 
momentum. Second, one can find a nonmonotonic increase in ((AiV) 2 ). Both 
these effects can easily be between one and two orders of magnitude greater 
than the statistical errors in present data. Third, one can find a nonmonotonic 
increase in the magnitude of (ApxAN). This quantity is small, and it has 
not yet been demonstrated that it can be measured. However, it may change 
at the critical point by a large factor, and is therefore worth measuring. In 
addition to effects on these and many other observables, it is perhaps most 
distinctive to measure the microscopic correlator (An p Ank) itself. The effects 
proportional to l/m^ in has a specific dependence on p and k. It introduces 
off-diagonal correlations in pk space. Like the off-diagonal anti-correlation 
introduced by energy conservation, this makes it easy to distinguish from the 
Bose enhancement effect, which is diagonal in pk. Near the critical point, the 
off-diagonal anti-correlation vanishes and the off-diagonal correlation due to 
sigma exchange grows. Furthermore, the effect of er exchange is not restricted 
to identical pions, and should be visible as correlations between the fluctuations 
of 7T + and 7r~. The dominant diagonal term proportional to 5 p k will be absent 
in the correlator (An+An^), and the effects of a exchange will be the dominant 
contribution to this quantity near the critical point. 

4 Pions From Sigma Decay 

Having analyzed the effects of the sigma field on the fluctuations of the direct 
pions, we next ask what becomes of the sigmas themselves. For choices of 
control parameters such that freeze-out occurs at or near the critical endpoint, 
the excitations of the sigma field, sigma (quasi)particles, are nearly massless 
at freeze-out and are therefore numerous. Because the pions are massive at 
the critical point, these er's cannot immediately decay into two pions. Instead, 
they persist as the temperature and density of the system further decrease. 
During the expansion, the in-medium a mass rises towards its vacuum value 
and eventually exceeds the two pion threshold. Thereafter, the it's decay, 
yielding a population of pions which do not get a chance to thcrmalize because 
they are produced after freeze-out. We estimate the momentum spectrum 
of these pions produced by delayed a decay. An event-by-event analysis is 
not required in order to see these pions. The excess multiplicity at low pr 
will appear and then disappear in the single particle inclusive distribution as 
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control parameters are varied such that the critical point is approached and 
then passed. 

In calculating the inclusive single-particle p-r-spectrum of the pions from 
sigma decay, we must treat m CT as time-dependent, and should also take G to 
evolve with time. However, the dominant time-dependent effect is the opening 
up of the phase space for the decay as m CT increases with time and crosses the 
two-pion threshold. We therefore treat G as a constant. We have estimated 
that in vacuum with m a = 600 MeV, the coupling is G ~ 1900 MeV, whereas 
at the critical end point with m a — 0, the coupling is reduced, perhaps by as 
much as a factor of six or so. In this section, we need to estimate G at the 
time when m a is at or just above twice the pion mass. We will use G ~ 1000 
MeV, recognizing that we may be off by as much as a factor of two. 

We parametrize the time dependence of the sigma mass by m a (t) = 
2m w (l + t/r), where we have defined t = to be the time at which m a has 
risen to Im^ and have introduced the timescale r over which m a increases 
from 2m w to Am^. It seems likely that 5 fm < r < 20 fm. We find that the 
mean transverse momentum of the pions produced by sigma decay is 

. . / 1000 MeV \ 2/3 /10 fm\ 1/3 
(pr)~0.58m^ J (_) . (8) 

We therefore estimate that if freeze-out occurs near the critical point, there 
will be a nonthermal population of pions with transverse momenta of order 
half the pion mass with a momentum distribution given in Ref. EJ. 

How many such pions can we expect? This is determined by the sigma 
mass at freeze-out. If m a is comparable to m„ at freeze-out, then there are 
half as many it's at freeze-out as there are charged pions. Since each sigma 
decays into two pions, and two thirds of those pions are charged, the result is 
that the number of charged pions produced by sigma decays after freeze-out is 
2/3 of the number of charged pions produced directly by the freeze-out of the 
thermal pion gas. Of course, if freeze-out occurs closer to the critical point at 
which m a can be as small as (6 fm) , there would be even more sigmas. We 
therefore suggest that as experimenters vary the collision energy, one way they 
can discover the critical point is to see the appearance and then disappearance 
of a population of pions with (j>t) ~ which are almost as numerous as 

the direct pions. Yet again, it is the nonmonotonicity of this signature as a 
function of control parameters which makes it distinctive. 

The event-by-event fluctuations of the multiplicity of these pions reflect 
the fluctuations of the sigma field whence they camel 3 ]. We estimate El that the 
event-by-event fluctuations of the multiplicity of the pions produced in sigma 
decay will be ((AN) 2 ) w 2.74(V). We have already seen in that the critical 
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fluctuations of the sigma field increase the fluctuations in the multiplicity of the 
direct pions sufficiently that the increase in the fluctuation of the multiplicity 
of all the pions will be increased by about 10 — 20%. We now see that in the 
vicinity of the critical point, there will be a further nonmonotonic rise in the 
fluctuations of the multiplicity of the population of pions with (px) ~ 
which are produced in sigma decay. 

5 Outlook 

Our understanding of the thermodynamics of QCD will be greatly enhanced 
by the detailed study of event-by-event fluctuations in heavy ion collisions. 
We have estimated the influence of a number of different physical effects, some 
special to the vicinity of the critical point but many not. The predictions of 
a simple resonance gas model, which does not include critical fluctuations, are 
to this point in very good agreement with the data. More detailed study, for 
example with varying cuts in addition to new observables, will help to fur- 
ther constrain the nonthcrmodynamic fluctuations, which are clearly small, 
and better understand the different thermodynamic effects. The signatures we 
analyze allow experiments to map out distinctive features of the QCD phase 
diagram. The striking example which we have considered in detail is the effect 
of a second order critical end point. The nonmonotonic appearance and then 
disappearance of any one of the signatures of the critical fluctuations which we 
have described would be strong evidence for the critical point. Furthermore, 
if a nonmonotonic variation is seen in several of these observables, then the 
maxima in all the signatures must occur simultaneously, at the same value 
of the control parameters. Simultaneous detection of the effects of the crit- 
ical fluctuations on different observables would turn strong evidence into an 
unambiguous discovery. 
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